subroutine comput_flux2
use var
implicit none

    integer m,j,i
    real(8) f(0:gridsx,0:gridsy,4)
    real(8) local_max
    real(8) f_m,f_p
    
    call comput_Qy_m
    call comput_Qy_p
    call L_F_local_y(Qy_m,Qy_p,f)
    do m=1,4
        do j=0,gridsy
            do i=1,gridsx
                flux2(i,j,m)=f(i,j,m)
            enddo                
        enddo  
    enddo
    
end subroutine